Edge State Magnetism of Single Layer Graphene Nanostructures 
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We study edge state magnetism in graphene nanostructures using a mean field theory of the 
Hubbard model. We investigate how the magnetism of the zigzag edges of graphene is affected by 
the presence of other types of terminating edges and defects. By a detailed study of both regular 
shapes, such as polygonal nanodots and nanoribbons, and irregular shapes, we conclude that the 
magnetism in zigzag edges is very robust. Our calculations show that the zigzag edges that are 
longer than three to four repeat units are always magnetic, irrespective of other edges, regular or 
irregular. We, therefore, clearly demonstrate that the edge irregularities and defects of the bounding 
edges of graphene nanostructures does not destroy the edge state magnetism. 
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I. INTRODUCTION 



Graphene is a single layer of carbon atoms forming a 
densely packed honeycomb lattice. Novoselov et al.— £ in- 
vented the top-down technique of isolating single layered 
graphene samples of a few microns in size by microme- 
chanical cleavage of graphite. This advance has lead to 
a flurry of activity, both theoretical and experimental, 
towards understanding the physics of graphene £is& 

A simple tight binding model^ for electron hopping on 
the honeycomb lattice produces two bands which touch 
each other at two points in the Brillouin zone. The chem- 
ical potential of the undoped graphene lies exactly at 
the energy where the two bands touch, implying that 
graphene is a zero gap system with two distinct "Fermi 
points" . What makes it even more interesting is the fact 
that the spectrum near these Fermi points resembles the 
Dirac spectrum of massless Fcrmionsj 7 -'2i 9 . and the density 
of states depends linearly on the energy. Graphene shows 
several interesting magnctotransport properties, such as 
quantum Hall effect at room temperature^ unconven- 
tional quantum Hall effect^ strongly suppressed weak 
localization magentoresistancei 2 - and quantum electrody- 
namics phenomena such as Klein paradox^ etc. Due to 
it's unusual electronic properties, graphene is a strong 
contender for future electronic applications. Examples 
are, graphene based field effect transistors (FET)^ spin 
valve devices,— gas sensors,— integrated ballistic carrier 
devices based on nanopatterned epitaxial graphene^ etc. 

Interesting nanostructures can also be made from 
graphene. For example, a nanoribbonil is obtained by 
reducing the dimension of a graphene sheet along one 
direction to the nanometric size, and a nanodot by re- 
ducing both the dimensions to the nanometric sizes. 
Graphene nanostructures can be terminated by many 
different types of edgesi 7 ., for example, by "zigzag" or 
"armchair" edges (see FIG. [TJ. These edges can have a 
profound influence on the electronic structure, and can 
give rise to interesting new phenomena. For example, 
zigzag edges have localized electronic states with nearly 
flat dispersion, giving rise to the finite DOS at the chem- 
ical potential, as has been reporte d 17 ' 18 from the the- 



oretical calculations. These have also been experimen- 
tally observed by scanning tunneling microscopy (STM) 
and spectroscop y 19 ' 20 . Indeed, the presence of the zigzag 
edges gives rise to unique physical and chemical proper- 
ties as reported by Son et al. — and Jiang et al. — based 
on the first principle density functional calculations. 

Focus of the present paper is on the edge state mag- 
netism of graphene nanostructures. Based on first princi- 
ple density functional calculations, it has been reported 
that the zigzag nanoribbons have an antifcrromagnctic 
(AFM) ground state, with one edge spin up and the other 
spin down 2 ^, whereas the armchair edged nanoribbons 
do not show any such magnetic property. The investi- 
gation has further been extended to graphene nanodots 
of rectangular shape 2 ^ and regular hexagonal, as well as 
triangular nanodots^ by first principle and mean field 
Hubbard model calculations; the results of these two ap- 
proaches are found to be in agreement £5. 

All the reports so far deal with the nanodots termi- 
nated by consecutive zigzag edges. However, the widely 
used top-down technique^ does not allow a careful control 
of the synthesis at the atomic level to produce graphene 
nanodots of regular shape, terminated, for example, en- 
tirely by zigzag edges. Furthermore, in spite of recent 
advances 2 ^, bottom-up techniques which are likely to 
give better control on the resulting structure are not yet 
widely practiced. The simplest defects to the "ideal" 
zigzag edged nanostructures are armchair terminations!^ 7 - 
Despite the fact that there are several studies of defect in- 
duced magnetism in graphene , 27 ' 28 ' 29 to the best of our 
knowledge, edge state magnetism of realistic graphene 
nanostructures, in particular those randomly terminated 
by the zigzag or armchair edges is yet to be examined in 
detail. In this work we address this issue and understand 
how the magnetism of zigzag edges "copes" with the pres- 
ence of other edges and random terminations. Using an 
unrestricted Hartree-Fock mean-field theory of the Hub- 
bard model on the graphene lattice, we study nanostruc- 
tures with various different terminations with the aim 
of uncovering how the edge magnetism is affected. Three 
cases are investigated. First, we study nanodots by start- 
ing from a perfect hexagonal shape and systematically 
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FIG. 1: (a) Hexagon: all the edges are of zigzag type. The size of the nanodot is denoted by N. We cut this polygon in 
different ways in order to produce various other polygons having both armchair and zigzag edges. Dark circles represent the 
edge sites in all the figures, and the dashed circles represent atoms that are removed, (b) Eight edged polygon: created by 
vertically cutting the hexagon resulting in two additional armchair edges per vertical side (quantified by N m ). The dotted 
portion in (b), (c) and (d) is the remaining part of the original full hexagon, (c) Ten edged polygon: obtained by cutting 
the hexagon along the slanted edge (no vertical cut) resulting in one armchair edge per slanted side (quantified by N n ) in the 
process, (d) Twelve edged polygon: created by simultaneously cutting along the vertical as well as the slanted edges, adding 
total four armchair edges of each type along the vertical and slanted directions respectively. 



introducing other types of edges. Second, we investigate 
graphenc nanoribbons terminated by the zigzag edges 
and systematically introduce "surface edge defects" and 
study the resulting magnetism. Third, we study "ran- 
dom" nanostructures by introducing completely irregu- 
lar edges, while preserving some zigzag segments at the 
edges. The Hartree-Fock approach allows us to study 
large nanostructures with as many as 1000 carbon atoms; 
such a study within the first principles framework would 
require much bigger computational resources. However, 
we do not take into account effects like bond length vari- 
ations and reconstructions at the edges, as well as effects 
of functional groups (e.g. oxidized edges) in the present 
calculation. 

The remainder of the paper is organized as follows. In 
the next section we brief the mean field Hubbard model 
used in our study. This is followed by a section (Sec. 
where we present and discuss our results. In the con- 
cluding section (Sec. HV)) . we discuss the significance of 
our results particularly in the context of experimental 
observation of edge magnetism in graphenc. 

II. MODEL DETAILS 

We study graphenc within a tight binding model in- 
cluding the on-site Hubbard repulsion, i. c., within the 



framework of the Hubbard mode l 25 ' 30 i 31 , which has been 
used earlier to study surface magnetisim 32 ' 33 The Hamil- 
tonian is 

h = -t ( c h c j° + h -°) + u Yl n it n n C 1 ) 

where Ci a (c\ a ) annihilates(creates) an electron of spin a 
at site i, rn a = c\ a Ci a is the number operator at site i, 
t is the hopping amplitude, U is the on-site Coulomb 
repulsion. The Hamiltonian has lattice translational and 
SU(2) spin rotation symmetries. Under the mean field 
approximation, we write the Hubbard Hamiltonian as 

H = -t^2 ( c h c j°+ h - c )+ u X^tKiH^i { n ir}-( n u}( n ir}) 

(ij)a i 

(2) 

where (rii a ) are the site occupations which are deter- 
mined self consistently. We perform a completely un- 
restricted calculation keeping all site occupations (n^) 
as unknowns. We focus on the undoped graphene, which 
had exactly one electron per site, i. e., we work with 
the half filled Hubbard model. From the self consis- 
tent ground state calculated numerically, wc obtain the 
expectation value of z-component of the spin operator 
Sf = hci a <Taa> c i<x' t° study the resulting magnetic struc- 
ture (a z is a Pauli matrix). 
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FIG. 2: (Online Color) Expectation value of (S z ) in four different graphene nanodots. The size of the arrows is proportional 
to the magnitude of magnetic moment. Note that there exist large magnetic moments at the zigzag edges compared to the 
sites lying inside or armchair edges. 



We first discuss known results 3 ^ of the magnetic phases 
of the half filled Hubbard model on an infinite honey- 
comb lattice. The honeycomb lattice does not break any 
symmetry of the Hamiltonian until the onsite Coulomb 
repulsion attains a critical valued of U c /t ~ 2.23. Above 
this value of the onsite Coulomb repulsion, the system 
breaks SU(2) spin rotation and lattice translational sym- 
metries and develops an antiferromagnetic order (due to 
its bipartite nature). This is unlike the 2D square lat- 
tice (with the nearest neighbour hopping), where even 
an infinitesimally small U will give rise to the antiferro- 
magnetic spin density wave(SDW)^ at half filling. This 
result of the honeycomb lattice may be understood start- 
ing from the magnetic susceptibility x(l) obtained from 
the linear response theory^ 



X(«) = (-1)^E 



f(e k -m) - /(£fe+ 9 - 1") 



£k — £*H 



(3) 



where /(e) is the Fermi function and q is the wave vector 
of perturbing field, is the electron dispersion, p is the 
chemical potential. Now the half filled Hubbard model 
on a bipartite lattice with nearest neighbour hopping has 
particle- hole symmetry. Using cqn. ([3]), the bare mag- 



netic susceptibility (for the antiferromagnetic response) 
can be obtained a a 35 ' 36 



(4) 



where p(e) is the density of states(DOS). In case of a 
square lattice, p{e) ~ log(e) in the vicinity of chemical 
potential at zero doping and xo has a logarithmic diver- 
gence. As a consequence, the generalized Stoner criterion 
for anti ferromagnetic (AF) susceptibility^ 



U > 



(5) 



is satisfied for even an infinitesimally small U and U c = 
for square lattice. However, in case of a honeycomb 
lattice, DOS depends linearly on the energy measured 
with respect to the chemical potential 3 ^, 



Pie) 



(6) 



Hence the divergence of Xo is absent (consult eqn. Q), 
i.e., xo is a finite quantity. Antifcrromagnctism, there- 
fore, appears only for Coulomb repulsion greater than a 
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critical value U c = 1/xo = 2.23i. Experimentally large 
graphene sheets do not show any magnetic ordering and 
hence the value of U in graphene should be less than 
U c (this result, of course, does not include fluctuational 
corrections). We shall use a value of U = 2t in our cal- 
culations, which at the mean field level will not produce 
any magnetic order in an infinite graphene sheet. 

As noted in the introductory section, zigzag terminated 
nanoribbons (quasi ID) possess special edge states which 
give rise to the flat bands and finite density of states at 
the Fermi level . 17 ' 18 ' 38 This leads to ferromagnetic spin 
orientations along a particular zigzag edge at infinitcs- 
imally small U c (from eqn. Q). On the other hand, 
two opposite zigzag edges, terminating the nanoribbon at 
two opposite sides, are found to be aligned antifcrromag- 
netically. The spin density (magnitude of the magnetic 
moment at a site) dies very quickly on moving into the 
"bulk" , i.e., the centre of the nanoribbon - this is indeed 
edge state magnetism. Armchair edges do not support 
any magnetic structure for values of U below U c , in con- 
trast to zigzag edges. 

In the next section we shall investigate various 
graphene nanostructures terminated by zigzag and other 
edges with the aim of investigating how the presence of 
other edges affects the edge state magnetism. 
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FIG. 4: (Color online) Average magnetic moment per edge 
site for the ten edged polygons. N„ is the number of armchair 
edges added along the slanted sides of the hexagonal polygon 
(FIG. [He)). Consult FIG. for the definition of <m z ), (m a ) 
and (m e ). Overall edge magnetization diminishes with the 
increase of number of armchairs. In particular, as we increase 
N n the horizontal zigzag segments are reduced in length and 
the most notable change occurs for N n greater than three, 
where the magnetization of the horizontal edges ({m z )) es- 
sentially vanishes. Note that for N n = 3 there are four repeat 
units in the horizontal zigzag segment and for N n = 4 the 
zigzag segments are two repeat units long. 



Before discussing each nanostructure in detail, we com- 
ment on the nature of the expected ground state. To- 
tal spin S z of the ground state of regular nanodots is 
found^ to be zero and this result is consistent with the 
predictions of Lieb's theorem^ This arises from the fact 
that the earlier calculations with regular nanodots will 
have equal number of atoms of the "A-sublattice" and the 
"B-sublattice" . On the other hand for irregular nanos- 
tructures, Na need not be equal to Nb, for example see 
FIG. [TJc) and (d) and total S z is expected (and found) 
to be nonzero in such cases. 



FIG. 3: (Online Color) Average magnetic moment per edge 
site for the eight edged polygons. N m is the number of arm- 
chair edges added along the vertical side (FIG. [TJb)). (m z ) 
is the average magnitude of the magnetic moment per site 
along the purely zigzag (horizontal sides in FIG. \Hh) and 
FIG.[2jb)) edges. (m a ) is the average magnitude of the mag- 
netic moment per site along the non-horizontal edges. (m e ) 
is the magnitude of the magnetic moment per site averaged 
over all the edge sites (zigzags as well as armchairs). As we 
increase N m , we add armchairs at the non-horizontal sides 
and hence (m a ), as well as (m e ), decreases. Moreover, in the 
process of increasing N m , we do not disturb the purely zigzag 
horizontal edges and (m z ) remains almost unaffected. 



III. GRAPHENE NANOSTRUCTURES: EDGE 
STATE MAGNETISM 

We study edge state magnetism in three types of 
nanostructures. Regular nanodots with zigzag and arm- 
chair edges, nanoribbons with "defected" zigzag edges, 
and random nanostructures with some zigzag segments. 



A. Nanodots 

Starting from a perfect hexagonal nanodot,— which 
is entirely enclosed by zigzag edges (FIG. Ufa)), we can 
make various different polygons, that are terminated by 
both zigzags and armchair edges (see FIG.QJb), (c),(d)). 
In this process we retain zigzag edges on the horizon- 
tal sides of the resulting polygon. Nomenclature of size 
etc. is explained in FIG. [TJ The sizes of the horizontal 
zigzag segments can by controlled by varying the lengths 
of the armchair segments N m and N n . 

The magnetic structure obtained from the mean field 
analysis of the Hubbard model is shown in FIG. [2] We 
observe that the surface sites, in particular those along 
the zigzag edges but not along the armchair edges, have 
large expectation values of the (S z ) operator. Moreover, 
the magnetic moment decreases sharply as we move to- 
wards the centre of the nanodot and at the "bulk" sites, 
it is in general at least one to two order of magnitude 
smaller in comparison to that at the zigzag edge sites. 
In all cases we find that the moment is along the "up 
direction" on sites of one sub-lattice and "down direc- 
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(b) 

FIG. 5: (Color online) Average magnetic moment per edge 
site for the twelve edged nanodots (FIG. Did)). See FIG. E] 
for the definition of (m z ) and (m a ). The size of the symbols 
is proportional to the magnitude of the average magnetic mo- 
ment, (a) The key finding is that the magnetic moment along 
the horizontal zigzag edges is affected by only their length 
and not by the presence of other edges. Increase of N m does 
not affect the length of the horizontal zigzag edges, while in- 
crease of N n decreases the length of the horizontal edges, (b) 
(m a ), average magnitude of magnetic moment along the non- 
horizontal sides decreases as we increase N m and N n , thereby 
adding more armchairs to the edges. 



tion" on the sites belonging to the other sub-lattice. We 
shall henceforth focus on the magnitude of the magnetic 
moment, keeping in mind this observation. 

We now present a quantitative analysis of the how the 
magnetism at the horizontal zigzag edges is affected by 
the introduction of other segments. FIG. O FIG. [Hand 
FIG. O show the plots of the magnitude of the magnetic 
moments per site for the three types of regular poly- 
gons constructed from a hexagonal nanodot discussed in 
FIG. [TJ The average magnitude of the magnetic moment 
per site along the zigzag segments is denoted by (m z ), 
while {m a ) denotes the average of the magnitude of the 
magnetic moment obtained on non-horizontal edges (the 
non-horizontal edges contain a mixture of both zigzag 
and armchair edges depending of N m and N n ). The 



average magnitude of the magnetic moment over all the 
edge sites is denoted by (m e ). FIG. [3] shows the varia- 
tion of various average magnetic moments for the eight 
edged nanostructure of shown in FIG.QJb). We find that 
the magnetization of the horizontal zigzag edges is unaf- 
fected by the presence of additional non-horizontal edges. 
The average magnitude of the magnetization (m e ) along 
all edges falls due to the presence of additional armchair 
segments. In FIG.[U we show the results of average mag- 
netization for the ten sided polygons. In this case the 
addition of non- horizontal armchair edges results (i. e., 
increase of N n ) in a reduction in the length of the hor- 
izontal zigzag segments. In this case we see that the 
magnetization of the zigzag segments falls (and so does 
the average magnetization) and for N n > 3, it essentially 
vanishes. The case of N n = 3 corresponds to having 
horizontal zigzag segments whose length is four repeat 
units. We turn now to the twelve sided polygons shown 
in FIG. [IJd) where we introduce two types nonhorizon- 
tal armchair segments denoted by N m and N n . As in the 
case of the eight sided polygon the length of the horizon- 
tal zigzag segment is unaffected by increase of N m while 
the length of the horizontal zigzag segments is reduced 
up on increase of N n . We see that the magnetization 
of the horizontal edges is unaffected by increase of N m , 
but strongly affected by increase of N n which reduces 
the length of the horizontal zigzag edges. The upshot of 
our calculation is that the zigzag edges retain a signif- 
icant amount of magnetization even in the presence of 
other edges as long as the length of the zigzag segments 
is greater than three to four repeat units. The number of 
repeat units of the zigzag edges that preserve the mag- 
netization depends on the value of U. For example, our 
calculations show that for U = 1.2t, the number of repeat 
units must be five to six for retaining the magnetization. 



B. Nanoribbons 

The nanodots discussed in the last section had "undis- 
turbed" zigzag edges. It is interesting to investigate if 
"defects" present in a zigzag edge can destroy the mag- 
netism. We investigate this issue by considering nanorib- 
bons which have a large aspect ratio. We start with 
nanoribbons enclosed solely by the zigzag edges, and in- 
troduce "armchair defects" along the long edges of the 
nanoribbon (see FIG. [5]). The number of armchair defects 
introduced is denoted by N a while the spacing between 
these defects is denoted by A^. Therefore if the value of 
Nd is smaller, the number of armchairs A^ is higher and 
vice versa. Here we define (m z ) as the average magnitude 
of the magnetic moment along the short edges, (m a ) as 
the average magnetic moment along the long edges possi- 
bly containing armchair defects. FIG. [6jd) shows a plots 
of these quantities. We observe from FIG. [H^d) that the 
moment along the short sides (which are defect free) is es- 
sentially unaffected by the number of defects on the long 
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FIG. 6: (Color online) Magnetic structure of graphene nanoribbons. Nanoribbon in (a) is enclosed by the zigzag edges. The 
long zigzag edges of the nanoribbons (b) and (c) is broken by intermittent armchair "defects". As we decrease Nd, the interval 
between the armchair "defects" (i.e., increase N a , the total number of armchair defects along the long side), the magnitude of the 
average localized magnetic moment along the long side (represented by (m a )) decreases. But (m z ), average localized moment 
along the short side remains unaffected. We denote average moment per edge site of the ideal zigzag structure (illustrated in 
(a)) by (m;) in (d). Note the sharp drop of (m a ) below Nd = 4, which is the limit of three zigzag repeat units. 



side. However (m a ) falls gradually with the decrease of 
Nd (increase of the number of the armchair defects) un- 
til Nd reaches about 4, where there is a dramatic fall 
in (m a ). Interestingly, this value of Nd corresponds the 
situation where the contiguous zigzag segments in the 
"defected zigzag edge" have about three to four repeat 
units. We therefore find a result that is fully consistent 
with that found in the case of regular polygon nanodots; 
as long as there are three to four repeat units of zigzag 
edge present, the ferromagnetism on the edge is robust. 



C. "Irregular" Nanostructures 



cation for the present study and focus on the nature of 
edge magnetism that arises in irregular structures. A 
study of FIG. [7] shows that the edge magnetism is very 
robust even in irregular structures; we find again that as 
long as there are three to four repeat units of a zigzag 
edge present, the edge magnetism is significant. We find 
an additional interesting feature. The sites with a single 
nearest neighbour (such as those that arc indicated by a 
solid arrow in FIG. [7]) have very large magnetic moments, 
larger than even the zigzag edges. As noted above, these 
are the "high energy edges" and are not likely to be sup- 
ported by thermodynamics. 



We shall now investigate several "irregular" nanodots 
created by "cutting out" patches out of a graphene sheet. 
Such a random cutout can produce edge atoms which 
have two nearest neighbours, or a single nearest neigh- 
bour. In the latter case the structure thus obtained may 
not be an equilibrium structure; we ignore this compli- 



IV. CONCLUSION 



As noted in the introduction there have been a large 
number of theoretical efforts investigating edge state 
magnctis m 23 i 24 i 25 i 40 i 41 , in regular nanostructures termi- 
nated by zigzag edges. However, we are not aware 
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FIG. 7: (Color online) Magnetic structure of irregular graphene nanodots. (a),(b) "Circular cut-out" nanodots. There are two 
different type of sites where the magnetic moment is high. First, the usual zigzag edges, lying along the horizontal side and 
marked by dashed arrow. Second, a set of sites, marked by solid arrows, have only one nearest neighbor, as opposed to two in 
case of zigzags and armchair edges. The second set of sites have even higher magnetic moment than compared to the zigzag 
sites. (c),(d) - nanodots with both regular edges and a "circular boundary". 



of any direct experimental observation of this phenom- 
ena for graphene nanoribbons or nanodots. There are 
a few experimental reports of magnetism in activated 
carbon fiber o 42 ' 43 , composed of a disordered network of 
nanographite, which remains to be the best proof of 
the theoretical prediction so far. In this paper we have 
enquired if "irregularity" of the experimental graphene 
nanostructures could be behind the lack of direct experi- 
mental corroboration. We have investigated how the edge 
state magnetism of the zigzag edges is affected by the 
presence of other edges, "defects" and random termina- 
tions. We find that the edge state magnetism is very 



robust to these "mutilations" of the nanostructures - for 
a value of U = 2t, as long as there are three to four re- 
peat units of a zigzag edge, the edge state magnetism 
is preserved. For smaller values of on site Coulomb re- 
pulsion, more number of repeat units of a zigzag edge 
is required for sustaining edge state magnetism (e.g., for 
U = 1.2t, the critical number of repeat units is five to 
six). In addition, we note that certain "high energy" 
edges (ones where the edge atoms have only one near- 
est neighbour) can have very large moments compared 
to even the zigzag edges. Thus our study demonstrates 
that the shape irregularity is unlikely to destroy the edge 
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state magnetism. However, it is clear that atomic res- 
olution magnetic force microscopy may have to be em- 
ployed for observation of edge state magnetism in short 
zigzag segment of graphene4^ On the theoretical front, 
the present approach is based on the mean field theory 
which, of course, neglects quantum fluctuation effects. 
These may be expected to be important in an "effec- 
tively ID" system like the zigzag edge of graphene and 
have to be investigated further. 

Note added: After completion of the manuscript, we 
became aware of a similar work by H. Kumazaki et al4£ 



We thank Kumazaki ct al. for bringing this to our notice. 
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